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The dynamics of an electron driven by arbitrary plane wave laser radiation is reformulated as a 
relativistically and mathematically exact refraction process based on an exact index of refraction. 
This reformulation leads to that index as an indicator of the energy transfer between the electron 
and the radiation. It also leads to the dynamics of the electron as being governed (i) by the 
Lorentzian version of what in Euclidean space is Snell's law, (ii) by an eikonal equation with the 
corresponding index of refraction, (iii) by geodesies on a spacetime manifold with in-general non- 
zero curvature ( "gemetrization of the laser radiation"), (iv) by the spacetime version of Fermat's 
principle of least time, and (v) by a Lorentzian stability criterion for the circumstance which in 
Euclidean space corresponds to the propagation of rays passing through a periodic wave guide of 
lenses which all have the same focal length. That criterion demands that the laser intensity satisfy 
laver < .56 ■ (10,000A/A)^ X 10^* watts/cm^ . 

PACS numbers: 41.75.Jv,45.,05.45.-a,45.05.-|-x,52.20.-j 



When a particle of charge q and mass m is placed into 
a laser beam whose radiation frequency uj /2t: has electric 
field amplitude Eq, then the particle executes oscillatory 
motion. The magnitude of the this effect is expressed by 
the dimcnsionless impulse factor 



qEo 



mcu! 



V ■ 



In light of (a) the simultaneous presence of an oscillat- 
ing magnetic field, and (b) the possibility of the motion 
being rclativistic, it is not surprising that the resulting 
complexity in the actual motion of the particle implies 
a corresponding complexity of the mathematical descrip- 
tion. 

However, it has turned out that, hidden behind this 
complexity, there often exists a readily identifiable sim- 
plicity, which physicists have expressed in terms of what 
is known as the "ponderomotive potential" and its gra- 
dient, the "ponderomotive force". These quantities arise 
from the fact that quite often the full motion of a par- 
ticle is characterized by two time scales. One charac- 
terizes the rapid quivering/oscillatory (fine-grained) as- 
pect of the full motion. The other characterizes a slower 
"guiding center" (coarse-grained) motion around which 
the particle executes its fast quivering oscillations [iHa] . 
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The ponderomotive force/potential, a slowly varying 
function of space and time, determines the slow motion 
dynamics of the particle. It is the result of performing 
a one-cycle average over the complex motion of the par- 
ticle. There are a number of ways[2-Q of doing this, 
but their common drawback is that one not only loses 
potentially useful information about the particle's mo- 
tion but, more importantly, misrepresents it, when the 
laser radiation becomes so intense {laver • (A/10, OOOA)^ > 
1.37 X 10^^ watts/cm <^ ?/ > 1) that the ponderomotive 
force varies as rapidly as the one due to the "rapid" quiv- 
erings/oscillations. Under such a circumstance, which 
includes a charge in a standing plane wavelJ,@,[3|: ^^ as- 
sumed decomposition into oscillatory plus averaged mo- 
tion along the direction of the laser beam, does not 
apply. (See, however, the second footnote below). For 
one thing, an ab initio averaging hides the possibility of 
resonance [3 where the frequency of the oscillatory mo- 
tion is a multiple of that of the averaged motion. 

Furthermore, femtosecond pulsed laser radiation 
makes any averaging scheme meaningless. There is not 
enough time to establish an average which changes slowly. 
Neither does averaging apply to charges located in, and 
hence scattered by, the transient overlap region of two 
counter propagating few cycle pulses. 

There is a superior way of understanding the dynamics 
of an electron driven by arbitrary plane wave laser radi- 
ation. We shall introduce an index of refraction which is 
exact for arbitrarily relativistic motion and/or arbitrarily 
short laser pulses and then reformulate the dynamics as 
a Lorentzian refraction process, as compared to its well 



known Euclidean cousin. Taking a cue from the prin- 
ciples of geometrical optics, this reformulation leads to 
(i) the Lorentzian version of what in Euclidean space is 
Snell's law, (ii) the spacetime version of Ferniat's princi- 
ple of least time, (iii) a Lorentzian stability criterion for 
the circumstance which in Euclidean space is the propa- 
gation of rays passing through a periodic wave guide of 
lenses which all have the same focal length. Although 
they are cousins, the Lorentz and the Euclidean refrac- 
tion processes are fundamentally distinct. Their respec- 
tive arenas, Euclidean space and space-time, entail (as 
shown below) opposite ray behavior, attraction for one, 
repulsion for the other. This difference can be traced to 
the difference in sign in the two respective Fermat varia- 
tional principles exhibited below. 

Dynamics of a Charge in a Generic Plane Wave Field 
- Consider an electron born at a generic location of the 
electromagnetic plane wave field of a laser directed along 
the z-direction. The four components of the vector po- 
tential are 

{Ao,Ai,A2,A3} = {0,0, A^{t,z),Ay{t,z)} . (1) 

The dynamics of the particle is governed by the 
four Lorentz (force law) equations of motion(9|, |lO| 
m(Px°'/dT^ = qF°l^dx^/dT {a = 0,1,2,3) for an elec- 
tron of mass m and charge q in the e.m. field F^^ = 
daAp — dfjAa. For a generic plane wave along the z- 
direction, these equations imply 






m 



<^{t,z) 



(fz d 

dr'^ dz 

dy 

m-— = Py - qAy(t,z) 

UT 

dx A r \ 

m—- =p^ - qA^(t,z). 
dr 



(2) 
(3) 
(4) 
(5) 



They are invariant under gauge changes having plane 
wave symmetry. Here time is measured in light traveling 
distance, {px,Py) is the particle's conserved ("canonical") 
transverse momentum, and 



Ht,^)^^U^-^A.At.. 



^-^A,(M) 
TO m 



(6) 
is the scalar potential which characterizes the longitudi- 
nal dynamics in the {t, z)-plane. 

The dynamics of the laser-accelerated charge is con- 
trolled entirely by Eqs.Q and ((3)). We therefore refer 
to them as the master system of equations. By contrast, 
Eqs. (Ul and ([5]) are merely slave equations. In relation 
to the longitudinal degrees of freedom these equations 
are dynamically passive. They have no effect on the par- 
ticle dynamics in the [t, z) -plane. Instead, they identify 
a physical measurable property, the transverse x and y- 
velocity components. 



The master system of equations implies that 



1 



1 

m 

2 



<^{t,z) = 'H^ 



(7) 



is an integral of motion. The integration constant ^ 
has been chosen such that the laboratory clock and a 
clock co-moving with the particle remain synchronized 

whenever the particle is at rest (^ = ^ = ^=0) in 

the lab frame. 

Dynamics Mathematized in Terms of Geometry - The 
two master equations lend themselves to being blended 
into a new mental unit, a geodesic on a Lorentzian man- 
ifold with a metric whose curvature is a geometrization 
of the electromagnetic field of the laser. This blending 
process is achieved in two steps: 

1. Note that these master equations. 



2^A 



d X 



„AB 



dT^ 



V 
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0; {2 



i:°^t,x^^z} (8) 



follow from the standard dynamical variational 
principle 



TO dx'^ dx'^ 
m dx dx^ 



lA, — ^dr 

TO dx" dx'' dx'' 



(9) 

dr. 
(10) 



A . 



Here upper case and lower case indices refer to the 
longitudinal {x 
coordinates {x" 
spectively. 



c" — t,x^ — z} and transverse 
y} coordinates re- 



^3 _ 



2. With the introduction of the transverse compo- 
nents Eq.Q-© and the addition of the inconse- 
quential term —pa^— the variational principle for 



the master equations becomes 



m dx dx^ . 



(11) 



In this variational principle with its two degrees of free- 
dom, the conserved transverse momentum components, 
{Pa '■ PxjPy}, are merely parameters, the transverse mo- 
mentum components of the particle. The dynamics of 
this system is governed by its Hamiltonian-Jacobi (H-J) 
equation 



1 



-V 



AB 



dS dS 



m 
$ = 

2 



2to, ' dx^ dx^ 
In geometrical optics the equivalent equation 
dS dS 



V 



AB 



dx^ dx^ 



2 2 

-m n 



(12) 



(13) 



with n^ given by Eg. ([251) below, is called the eikonal equa- 
tion llf having index of refraction n. 



The dynamics of our system on the 2-d (flat) 
Minkowski spacetime is a study of those particles hav- 
ing rest mass which varies with spacetime location as 
specified by the e.m. field on the right side of the H- 
J equation. However, such a study is mathematically 
equivalent to the one of particles with fixed rest mass 
but on a manifold with a metric which is location depen- 
dent. This statement is expressed mathematically by the 
H-J equation 



AB dS dS 2 n 



dx^ dx^ 



whose inverse metric 



.AB 



^AB 






(14) 



(15) 



is location dependent, and whose corresponding metric 
is HP 



9AB 



1- 

2(7 



\m m / \m m 



VAB 
= Vabs 
Its curvature is [IJ, ll3| 

R^BCD^T^iScVBD 



5^ 



dVbc) 



(16) 
(17) 

(18) 



This means, among other things, that the curvature van- 
ishes whenever the laser radiation {A^, Ay} is such as to 
satisfy 



n. 






dz 



= 



which is to say, the radiation propagates into the +z or 
the —z direction, while any mixture of the two will result 
in non-zero curvature. In the first case the system is 
integrable; the H-J equation can be solved by the method 
of separation of variables. In the second case the laser- 
particle system is one which is non-integrable. 

Laser Driven Dynamics in a Longitudinal Electric 
Field A grasp of laser-induced charge dynamics on a 
fundamental level requires taking cognizance of strong 
Coulomb ("space charge") fields. Although their forces 
are strictly longitudinal, their indirect cause is the strong 
charge separation due to the ultra-intense laser interac- 
tion with the target environment |14| - [l7| . 

With such fields present, instead of Eq.(IT|), the com- 
ponents of the vector potential are 



d'^x'^ 



-,AB 



a$ 



V 



= qv 



„AB 



F, 



dx'' 



BC- 



(21) 



dT^ ' ' dxB "■' ^ ^" dT 

The reduced variational principle giving rise to them is 



m dx dx^ dx^ 



$ Wr. 



(22) 



and the corresponding geomctrized H-J equation is 



AB f dS , \ f dS 



m^ = 



(23) 



Its inverse metric is still given by Ea. ([T5)) . 

Dynamics as Geometrical Optics - The integral of mo- 
tion, Eq. ([7]), has a property that leads directly to a 
geometrical optics formulation of the dynamics of the 
particle: Rewrite the integral in the form 



= n'^{t,z) 



(24) 



where 
n(t, z) 



P^-±A4t,z)) +(^-AA,(t,z; 

mm / \m m 



(25) 



The integral, Eq. ((24| . holds along every spacetime tra- 
jectory. This suggests that one should replace the proper 
time T of the charge with its longitudinal proper timc^ 



df 



n (^(t), z{t)) dr 



(26) 



as the world line parameter. Such a replacement yields 
^ ~ n (> 0). Consequently, the integral of motion be- 
comes 



1 



(27) 



and Eqs. ([2])- ([3]) become the optical master equations 
d dt 



dr dr 
d dz 

dr df dz 



dn 

"dt 
dn 



(28) 
(29) 



or more succinctly. 



{Ao,Ai,A2,A3} = {At(t,z),A,{t,z),A^{t,z),Ay{t,z)} , 

(19) 

while the master Eqs.® with its additional Coulomb 
field 



d dx^ 
df df 



AB 



dn 



^AB 



1 
-1 



A 
B 



= 0,1 
(30) 



Fbc = dBAc - dcAs 



(20) 



^ This time is measured by a clock moving along fixed x and y 
coordinates, but co-moving (and accelerating) with the charge 
strictly along the z-direction. 



Compare these equations with those for a hght ray 
through a medium with refractive index n(x^,x^^x'^) in 
geometrical optics p^. 



d dx^ _ „j dn ^ i 
ds ds dx^ ' J 



1,2,3, 



(31) 



where i5*^ is the Kronecker delta, a sum over repeated 
indices is understood, and s is the geodesic length pa- 
rameter along the ray. One sees that if n{x^,x'^,x^) is 
the familiar Euclidean index of refraction for a medium 
in Euclidean space, then Eq. (j25p is the Lorentzian index 
of refraction for the e.m. -induced medium in Lorentzian 
spacetime. The former determines the ray trajectories 
a;'(s); the latter determines the world lines x^{f). 




- z-displacement along laser beam - 



FIG. 1: World line of a charge through the periodic spacetime 
lattice of the electromagnetic standing wave field of a laser. 
The periodic index of refraction, Ea. (|35|) . has maxima which 
are surrounded by the oval shaped isograms. The gradient 
vector field is (Lorentz) orthogonal to these isograms. Fur- 
thermore, this vector field pushes the charge and thereby im- 
parts to it the acceleration given by Eq. (|28|I - P9|I . For 77 = 1.5 
and for the depicted initial conditions a numerical integra- 
tion of these equations yields the world line reproduced in 
this figure. 

The optical master equations are a local expression of 
5 j n{t, z)\/ dt"^ — dz^ = , which is what in Euclidean 
space corresponds to Format's principle of least time, 
6jn{x\x^,x^)^{dx^)^ + {dx^y + {dx^Y = . 

All world lines are determined by the optical master 
Ea. ((28|) - ((29)) . A cursory examination reveals that Eq. ((27)) 
together with Eq.([28l) (resp. Eq.(l29l) ) implies Eq.(|29|) 
(resp. Eq. (|28| ). Consequently, one readily finds that the 
world lines satisfy the following single differential equa- 
tion. 



a z 
'dp 



1 - 



1 
2^ 



dn'^ dz dn^ 
dt dt dz 



(32) 



Given any laser-induced refractive index, Eq. (j25p . this 
master equation governs all properties of all possible elec- 
tron trajectories resulting from a generic plane wave laser 
environment. 



One of the properties of a spacetime trajectory is how 
it gets bent (a.k.a. accelerated) by the refractive index 
gradient. The right hand side of Eq.([32|) is proportional 
to the directional derivative along the co-moving z-axis. 
Indeed, given that, in light of Eq. ((27)) . the unit tangent 
to the 2-d world line is (^j^), the unit vector along 
the co-moving z-axis is (||, ^) = z. It follows that 



' dz\ 
df \dtj 



dz\ ' 
dt) 



dz Inn 



(33) 



where dzlnn = ^dtlnn + -^dzlnn. Thus, if n in- 
creases along the positive (resp. negative) z-direction, 
the charge will accelerate into the negative (resp. pos- 
itive) z-direction. In other words, opposite to the Eu- 
clidean circumstance, regions of tiigher Lorentzian index 
of refraction accelerate charges into a direction where the 
index is smaller. A charge gets expelled. 

One manifests an instance of particular interest by 
stipulating the vanishing of the particle's conserved 
transverse momentum: {px,Py) = (0,0). Under this 
stipulation, Eq. ipS)) indicates that the square of the 
Lorentzian index is proportional to the field intensity. 
Hence one may conclude that charges drift away from re- 
gions of locally higher intensity regardless of the signs of 
the charges, a well-known phenomenon which is usually 
attributed to the presence of the so-called ponderomotive 
force. 

A second property arises for a particularly interest- 
ing case, a charge driven by the laser field of a standing 
wave(j,|6|,LZ| linearly polarized along, say, the x-axis. For 
such a wave the e.m. vector potential is 

Ax{t,z) = — sinujtsmLuz , Ay{t,z) =0 . (34) 

and the squared index of refraction is 

n {t, z) = 1 + {px/rn — i]iimujtsmujz) . (35) 

Such a refractive index n{t, z) forms a two-dimensional 
optical lattice of what in Euclidean space are convex and 
and concave lenses. Referring to Figure'^ [U one infers 
from the depicted trajectory that there exists a linear 
stack of of convex "Lorentzian" lenses all localized around 
z = half way between two adjacent maxima of the re- 
fractive index n(t,z). In Figure [T] they are separated 
("vertically") by an amount L = tt. These lenses guide 
the center of the oscillating trajectory. Moreover, each 
lender has a temporal focal length. For a laser beam hav- 
ing a standing wave field of frequency j^ , this focal length 
is 



nu} rf 



(36) 



^ In this figure we have set px = 0. This suppresses certain in- 
teresting features associated with averaged motion, which have 



been identified in[y,[3|. 



It is a property pertaining only to "paraxial" trajectories, 
i.e. those for which the lab velocity of the charge is non- 

relativistic: (^) < 1. 

The mathematical reasoning that leads to Ea.lfM]) is 
taken directly from geometrical optics. It starts with 
Eq. (f32|) combined with Eg. ([35]) . The idea is to calculate 
for a single lens 



-i§ 



I-KJU 



dt^ 



dt 



the change in particle velocity over time 27r/w for trajec- 



tories near z = that start with 
such paraxial spacetime trajectories is 



dz 
dt 



0. The result for 



dt F 



where F^ the temporal focal length, is given by Eg . (p6| . 
Stable and Unstable Motion - The existence of a focal 
length and the periodic structure encountered by a parax- 
ial particle trajectory direct attention to the possibility 
of parametric instability in its oscillatory motion. Recall 
that from geometrical optics one knows that if neighbor- 
ing convex lenses of a periodic stack are separated by the 
same amount L, then such a lens system accommodates 
linearly stable paraxial trajectories if and only if [13, 1201 

[Separation length) _ L 
[Focal length) F 

From Figure [1] one sees that the separation between con- 
secutive lenses is 

c 

u 

In light of Eq. ((36l) one finds that paraxial spacetime tra- 
jectories oscillate around z = Q, i.e. are linearly stable if 
and only if the laser impulse parameter satisfies 






(37) 



which is to say, the laser intensity should satisfy laver < 
.56- (10,000A/A)2 X lO^Vatts/cm^. If this inequality is 
violated one has a possible type of parametric resonance, 
which we alluded near the beginning of this article. For a 
differential geometric and more detailed analysis of how 
such a parametric resonance may lead to the breakdown 
of the ponderomotive approach, one is referred to [l^] • 

Snell's _Lau;[l3| - A fundamental manifestation of a 
refraction process is the manner in which a ray propa- 
gates across the boundary between two regions having 
different indices of refraction. In Euclidean space this 
propagation is expressed by Snell's law. What is its form 
for a spacetime medium whose index of refraction varies 
(in the limit) discontinuously, as in Figure [2F 

Consider the world line of a charge as it leaves a space- 
time medium with index n_ and enters another one with 



index n_|_ . Fig. [2] gives a close-up view. One introduces 
the "null" coordinates 

u ^ t — z "retarded time" coordinate (38) 

V ^ t + z "advanced time" coordinate (39) 

and lets the boundary be the locus of events where 

V = "history of pulse discontinuity" . (40) 

It is along characteristics like this where solutions to the 
wave equation 



d^ f A^{u,v) 



dudv \Ay{u,v) 



(41) 



are allowed to be discontinuous, and hence where the 
index of refraction, Eq. (PS|) . is allowed to be discontin- 
uous. 

A charge which crosses such a boundary will experience 
a refractive index of the form 



n{u, v) 



?7_ u < 
n+ < f 



(42) 



The indices are different but constant on either side of 
the boundary. 




(space) 



FIG. 2: Refraction of the world line of a charge 

The problem is indicated in Figure [51 One must estab- 
lish the relationship between the slopes 



dv _d{t + z) _l + dz/dt 
du d{t — z) 1 — dz/dt 



(43) 



on either side of the "null" boundary ( "characteristic" of 
the wave equation) v = 0. The solution consists of the 
statement that 



dv 
du 



dv 
du 



(44) 



where the "4-" and "-" refer to < w and v < respec- 
tively. This is the Lorentzian version of what is Snell's 
law in Euclidean space. 

The slope dv/du is well known. It is the square of the 
Doppler factor ("rapidity factor") e^ for a particle with 
z-velocity ^ = tanhe* = (e^-e"^)/(e^-|-e"^). The quan- 
tity 9 is generally known as the particle's rapidity. The 
differences and similarities between the Lorentzian and 
the Euclidean versions of Snell's law become particularly 



perspicuous when one introduces this Doppler factor. In 
terms of it Snell's law takes the formfTj 



k'>0 



e n 



v<0 



(45) 



across a left-travehng pulse (v = 0). This is to be con- 
trasted with the Euclidean version of Snell's law, which, 
we recall, is 



ni sm f/i = 71,2 sm 6^2 



(46) 



The validity of Eq. (|44]), and hence of Eq. (gS]), follows 
directly from Ea. (|27|) and the optical master Eqs. ((28)) - 
(l29l) recast in terms of the "null" coordinates u and v. 
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Fij 



= 2 



dn{u, v) 
du 



(47) 



(48) 



f with II as a worldline parameter. With the help of the 
second equation of Eq. (|48| . Eq. (jT7|) leads directly to 
Eq. dH. 

Refractive Index as a Measure of Kinetic Energy Gain 
- Recall that the charge kinetic energy is given by 



K{t) 



dt 
dr 



Thus, along a given trajectory, 



dK 
d7 



m d 

Tat 



n^{t,z) 



(49) 



(50) 



by Eqs. ©, ©, and (P5|) . This equation implies that a 
charge gains kinetic energy as it enters a region of higher 
refractive index from one where the index is smaller. 
Conversely, the charge loses kinetic energy as it leaves 
a region of higher refractive index and enters one where 
the index is smaller. In other words, a change in the 
Lorentzian refractive index is an indicator of energy ex- 
change between the charge and the e.m. field. 

Furthermore, since regions of higher index accelerate 
charges into a direction where the index is smaller, it fol- 
lows that charges are inclined to lose their kinetic energy 
in response to the e.m. field. 

The relation between the kinetic energy gain and the 
refractive index is most transparent for charges respond- 
ing to a traveling wave whose vector potential depends 
only on, say, the advanced time coordinate u = t — z. In 
this case, Eq. ([50)1 becomes 



dK 

du 



2pu du 



(51) 



where p^ is the conserved value of du/dr. Consequently, 
the kinetic energy gain changes linearly with the increase 
in the square of the refractive index. That is. 



AK 



TTiAn^ 
2pu 



(52) 



Energy Gain from. Interaction with a Chirped Traveling 
Pulse ~ As an application, consider the particle dynamics 
in the chirped plane e.m. field which is polarized along 
the a;-axis {Ay = 0) and whose field function is given by 



dA^ 
du 



-£ cos ((/?o + u + bu^) g{u) 



(53) 



as in [2l|, |22|. Here £, ipo, b, and u are the field ampli- 
tude, a constant initial phase, and a dimensionless chirp 
parameter, and the dimensionless retarded time coordi- 
nate u = wqu, where wq is the frequency of the field at 
(t,z) = (0,0). Also, 



liu) 



exp 



(u - 4a) 



21 



2a2 



(54) 



is a pulse-shape function, where a is related to the 
pulse duration r (full-width at half-maximum) via a = 

one deduces immediately the 



(55) 



woT/(2\/21n2). From §^ 

evolution of the particle kinetic energy in u: 



AK = 



2mpu 



cos {ipQ + u + bu ) g{u) du 



which matches the calculation in [2l|, |22|. Specifically, 
one sees that the particle energy gain scales linearly with 
the field intensity as a consequence of the linear relation 
between the energy gain and the change in the refractive 
index squared. 

Conclusion - With a Lorentzian index of refraction at 
its foundation, relativistic particle dynamics driven by 
generic plane wave laser radiation of arbitrary intensity 
has been mathematized in terms 

1. of the Lorentzian ( "spacetime" ) version of the op- 
tical master equations, 

2. of geodesies on a spacetime manifold with in- 
general non-zero curvature ( "geometrization of 
laser radiation"), 

3. of Fermat's principle of least time, 

4. of Sneh's law, 

5. of the focal lengths of Lorentzian lenses, which 
make up the periodic spacetime lattice of a stand- 
ing wave, and 

6. of the associated stability criterion for particles 
moving through this lattice, 

the first four of which are mathematically equivalent. 
Last, but not least, the Lorentzian index of refraction is 
shown to be an indicator of energy transfer between the 
radiation and the particle. 



[1] 

[2] 
[3] 
[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 



L. Landau and E. Lifschitz, Mechanics (Pergamon Press 
Inc., Elmsford, NY, 1969), pp. 93-95, 2nd ed. 
T. Kibble, Phys. Rev. Lett. 16, 1054 (1966), U RL 
http://link.a ps . org/doi/10 . 1103/PhysRevLet t . 16. 1054[ 
T. Kibble, P hys. Rev. 150, 1060 (1966), URL 
|http : //link. aps . org/doi/10 . 1103/PhysRev . 150 . 1060| 
D. Bauer, P. Mulser, and W. H. Steeb, 
Phys. Rev. Lett. 75, 4622 (1995), URL 
|http : //link . aps . org/doi/10 ■ 1103/PhysRevLett . 75 . 4622 [ 



[12] 
[13] 



E. A. Startsev and C. 
Phys. Rev. E 55, 7527 



J. McKinstrie, 
(1997), URL 



|http : //link . aps . org/doi/10 . 1103/PhysRevE . 55 . 7527 [ 
A. E. Kaplan and A. L. Pokrovsky, 
Phys. Rev. Lett. 95, 053601 (2005), URL 



[14] 
[15] 
[16] 

[17] 



[18] 



|http : //li nk . aps . org/doi/10 . 1103/PhysRevLett . 95 . 05360 
A. L. Pokrovsky and A. E. Kaplan, [19] 
Phys. Rev. A 72, 043401 (2005), URL 
http://link.ap s . org/doi/10 . 1103/PhysRevA . 7 2.043401] [20] 
V. L Arnold and A. Avez, Ergodic Problems of Classical 
Mechanics (W.A. Benjamin, Inc., New York, NY, 1968), [21] 
pp. 81-109. 

C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravita- 
tion (W.H. Freeman and Company, San Francisco, 1973), 
chap. 3.1. [22] 

J. Jackson, Classical Electrodynamics (John Wiley and 
Sons, Inc., New York, NY, 1975), pp. 551,572, 2nd ed. 
M. Born and E. Wolf, Principles of Optics (Pergamon 



Press Inc., Elmsford, NY, 1980), p. 112, sixth ed. 
K. H. Kwa, J. Phys. A-Math. Theor. 45, 105203 (2012). 
K. H. Kwa, Ph.D. thesis, Dept. of Math- 
ematics, O.S.U., Dept. of Mathemat- 
ics, O.S.U., Columbus, OH (2009), URL 
Ihttp : //rave . ohiolink . edu/etdc/view? acc_num=osul250 103994 [ 
F. Brunei, Phys. Rev. Lett. 59, 52 (1987). 
F. Brunei, Phys. Fluids 31, 2714 (1988). 
J. P. Geindre, R. S. Marjoribanks, and P. Audebert, 
Phys. Rev. Lett. 104, 135001 (2010). 
M. Veltcheva, A. Borot, C. Thaury, A. Malvache, 
E. Lefebvre, A. Flacco, R. Lopez-Martens, and V. Malka, 
Phys. Rev. Lett. 108, 075004 (2012). 
M. Born and E. Wolf, Principles of Optics (Pergamon 
Press, Inc., Elmsford, NY, 1980), chap. 3.2, 6th ed. 
D. Marcuse, Light Transmission Optics (Van Norstrand 
Reinhold Inc., New York, 1972), pp. 174-184. 

A. Yariv, Quantum Electronics (John Wiley and Sons, 
Inc., New York, 1975), pp. 99-104, 2nd ed. 

B. J. Galow, Y. I. Salamin, T. V. Liseyk- 
ina, Z. Harman, and C. H. Keitel, Phys. 
Rev. Lett. 107, 185002 (2011) , URL 
|http : //link . aps . org/doi/10 . 1103/PhysRevLett . 107 ■ 185002[ 
Y. I. Salamin, Phys. Lett. A 376, 2442 (2012), URL 
|http : //dx ■ doi . org/10 ■ 1016/ j ■ physleta . 2012 . 06 ■ 020 [ 



